EXISTENCE ANALYSIS OF MAXWELL- STEFAN SYSTEMS FOR 
MULTICOMPONENT MIXTURES 



ANSGAR JUNGEL AND INES VIKTORIA STELZER 

Abstract. Maxwell-Stefan systems describing the dynamics of the molar concentrations 
of a gas mixture with an arbitrary number of components are analyzed in a bounded 
domain under isobaric, isothermal conditions. The systems consist of mass balance equa- 
tions and equations for the chemical potentials, depending on the relative velocities, sup- 
plemented with initial and homogeneous Neumann boundary conditions. Global-in-time 
existence of bounded weak solutions to the quasilinear parabolic system and their expo- 
nential decay to the homogeneous steady state are proved. The mathematical difficulties 
are due to the singular Maxwell- Stefan diffusion matrix, the cross-diffusion coupling, and 
the lack of standard maximum principles. Key ideas of the proofs are the Perron- Frobenius 
theory for quasi-positive matrices, entropy-dissipation methods, and a new entropy vari- 
able formulation allowing for the proof of nonnegative lower and upper bounds for the 
concentrations. 



1. Introduction 

The Maxwell-Stefan equations describe the diffusive transport of multicomponent mix- 
tures [221 121] • Applications include various fields like sedimentation, dialysis, electrolysis, 
ion exchange, ultrafiltration, and respiratory airways [El [27]. The model bases upon inter- 
species force balances, relating the velocities of the species of the mixture. It is well-known 
that the usual Fickian diffusion model, which states that the flux of a chemical substance 
is proportional to its concentration gradient, is not able to describe, e.g., uphill or os- 
motic diffusion phenomena in multicomponent mixtures, as demonstrated experimentally 
by Duncan and Toor [13]. These phenomena can be modeled by using the theory of non- 
equilibrium thermodynamics, in which the fluxes are assumed to be linear combinations 
of the thermodynamic forces [HI Chap. 4]. However, this model requires the knowledge of 
all binary diffusion coefficients, which are not always easy to determine, and the positive 
semi-definiteness of the diffusion matrix. The advantage of the Maxwell- Stefan approach is 
that it is capable to describe uphill diffusion effects without assuming particular properties 
on the diffusivities (besides symmetry). 
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We consider an ideal gaseous mixture consisting of iV + 1 components with molar con- 
centrations Ci(x, t) for i — 1, . . . , N + 1 (see Appendix |AJ. Since we concentrate our study 
on cross-diffusion effects, we suppose isothermal and isobaric conditions. Then the total 
molar concentration Xlilt 1 Ci is constant and we set this constant equal to one. More gen- 
eral situations are investigated, e.g., in [15J. The dynamics of the mixture is given by the 
mass balance equations 

(1) d t Ci + div Ji = fi(c) in Q, t > 0, % — 1, . . . , N + 1, 

where Ji = CiU% denotes the molar flux, Ui the mean velocity, the net production rate of 
the 2-th component, c = (c±, . . . , cn+i) t , and Q C M d (d > 1) is a bounded domain. We 
have assumed above that the averaged mean velocity vanishes, J^^ 1 CjUj = 0. Then the 
conservation of the total mass implies that the total production rate vanishes, X^t 1 r «( c ) = 
0. The fluxes are related to the concentration gradients by 

N+1 r T -r T 

(2) Vc i = - D . > i = l,...,tf+l, 

where Dij > for i ^ j are some diffusion coefficients. The derivation of these relations is 
sketched in Appendix |A] 

The aim of this paper is to prove the global-in-time existence of weak solutions to 
system (DO)-© for constant coefficients > 0, supplemented with the boundary and 
initial conditions 

(3) Vci-z/ = ondQ, t > 0, c 4 (-,0) = c° in Q, i = 1, . . . , iV + 1, 

where v is the exterior unit normal vector on dQ. There are several difficulties to overcome 
in the analysis of the Maxwell-Stefan system. 

First, the molar fluxes are not defined a priori as a linear combination of the concen- 
tration gradients, which makes it necessary to invert the flux-gradient relations (F2J). As 
the Maxwell-Stefan equations are linearly dependent, we need to invert the system on a 
subspace, yielding the diffusion matrix A' 1 . In the engineering literature, this inversion is 
often done in an approximate way. For instance, a numerical solution procedure for N = 3 
was developed in [2] and the special case = for some constants fi > was 

investigated in [3]. 

Second, equations ([I])-© are coupled, which translates into the fact that A~ l is generally 
a full matrix with nonlinear solution-dependent coefficients. Thus, standard tools like the 
maximum principle or regularity theory are not available. In particular, it is not clear how 
to prove nonnegative lower and upper bounds for the concentrations. Moreover, it is not 
clear whether A' 1 is positive semi-definite or not, such that even the proof of local-in-time 
existence of solutions is nontrivial. 

Third, it is not standard to find suitable a priori estimates which allow us to conclude 
the global-in-time existence of solutions. 

In view of these difficulties, it is not surprising that there are only very few analytical 
results in the mathematical literature for Maxwell- Stefan systems. Under some general 
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assumptions on the nonlinearities, Giovangigli proved that there exists a unique global 
solution to the whole-space Maxwell-Stefan system if the initial datum is sufficiently close 
to the equilibrium state [I5j Theorem 9.4.1]. Bothe |1] showed the existence of a unique 
local solution for general initial data. Boudin et al. [7] considered a ternary system (N = 2) 
and assumed that two diffusivities are equal. In this situation, the Maxwell- Stefan system 
reduces to a heat equation for the first component and a drift-diffusion-type equation for 
the second species. Boudin et al. [7J proved the existence of a unique global solution and 
investigated its long-time decay to the stationary state. Up to now, there does not exist a 
global existence theory for fll}-© for general initial data. We provide such a result in this 
work. 

After inverting the flux-gradient relations (J2J) on a suitable subspace, the Maxwell- Stefan 
equations become a parabolic cross-diffusion system. Amann derived in [TJ sufficient condi- 
tions for the solutions to such systems to exist globally in time. The question if a given local 
solution exists globally is reduced to the problem of finding a priori estimates in suitable 
Sobolev spaces. Another approach was developed in [SJ [HI (TTJ [19] for systems arising 
in granular material modeling, population dynamics, cell biology, and thermodynamics to 
treat cross-diffusion systems whose diffusion matrix may be neither symmetric nor positive 
semi-definite. The idea is to exploit the entropy structure of the model by introducing 
so-called entropy variables. In these variables, the new diffusion matrix becomes positive 
semi-definite and an entropy-dissipation relation can be derived. However, in all of the 
mentioned papers (except [9]) systems with two equations only have been considered. 

In this paper, we combine and extend the entropy-dissipation technique of [HI HI HU (TTJ 
IT9] as well as ideas of Bothe [I] to overcome the above mathematical difficulties. We are 
able to prove the global-in-time existence of weak solutions to ([I])-© for arbitrary diffusion 
matrices and general initial data. This result is obtained under the following assumptions: 

• Domain: Slc» d (d<3) is a bounded domain with d£l G C ' . 

• Initial data: c°, . . . , c% (N > 2) are nonnegative measurable functions, c° N+1 = 

N 

(4) E^ 1 - 

i=l 

• Diffusion matrix: (-Dy) € ]R( Ar+1 ) x ( Ar+1 ) is a symmetric matrix with elements Dy > 
for i 7^ j. 

• Production rates: The functions r.; G C°([0, l]^ 4 " 1 ; K), i — 1, . . . , N + 1, satisfy 

N+l N+l 

(5) yVi(c) = 0, ^ ri(c) log q <0 for all < ci, . . . , cjv+i < 1- 

i=l i=l 

We stress the fact that, although the diffusion coefficients Dij are constant, the diffusion 
matrix of the inverted Maxwell-Stefan system depends on the molar concentrations in a 
nonlinear way (see below) and we need to deal with a fully coupled nonlinear parabolic 
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system. Our proof also works for diffusion coefficients depending on the concentrations c« 
if the coefficients d^ are bounded from above and below. 

The regularity on the boundary dQ is needed for the a priori estimate ||io||.ff 2 (n) < 
C||/l|L 2 (n) °f the elliptic problem —Aw + w — f in ft, Vw ■ v = on dQ. 

The inequality imposed on the production rates is needed to prove that the entropy is 
nonincreasing in time. It is satisfied if, for instance, N = 4 and r% = r% = C2C4 — C1C3, 
r 2 = r 4 = C1C3 — C2C4 [10]. For the existence result, the inequality can be weakened by 

N+l 

(6) y"Vj(c) logQ < C r for all < ci, . . . ,cjv+i < 1, 

i=l 

where C r > is some constant independent of Cj (see ff27|) ). This condition is satisfied, for 
instance, for the tumor-growth model in [T9] . 

Our first main result is the global existence of solutions to (TTJ)-([3]). 

Theorem 1 (Global existence of solutions). Let the above assumptions hold. Then there 
exists a weak solution (c±, . . . , cn+i) to ([I])-© satisfying 

Ci G Lf oc (0, 00; H\Q)), dta G LL(0, 00; V), 

N 

< Ci < 1, i = 1, . . . , N, cn+i = 1 — Ci > in O, t > 0, 

i=l 

where V is the dual space ofV— {u G H 2 (Q) : Wu ■ v = on dQ}. 

To be precise, the existence theorem for ([I])-© has to be understood as an existence 
result for a system in N components (see below), which is equivalent to (JT|)- (J3j) as long as 
q > for alH = 1, . . . , N and J2f=i c « < 1 are satisfied. 

We explain the key ideas of the proof. For this, we write (|2j) more compactly as 

(7) Vc = A(c)J, 

where A(c) G R(N+iMN+i) and Vc = (da/dx^n, J = (J u . . . , Jjv+i) t G R^ N+1 ^ d . Using 
the Perron-Frobenius theory for quasi-positive matrices, Bothe [1] characterized the spec- 
trum of A(c) in case that Cj > for i — 1, . . . , N and J2f=i c i < 1- Under these conditions, 
A(c) can be inverted on its image. Then, denoting its inverse by A(c) -1 , ([Tl)-([2|) can be 
formulated as 

(8) d t c - div(-i(c)~ 1 Vc) = r(c) in fi, t> 0, 

where r(c) = (ri(c), . . . , rAr +1 (c)) T . 

It turns out that it is more convenient to eliminate the last equation for cjv+i, which is 
determined by Cat + i = 1 — J2iLi c h anc ^ to work only with the system in N components. 
We set d = (c u . . . , c N ) T , J' = ( J u . . . , J N ) T , and r'(c) = (r x (c), . . . , r^c))^ Using 

the facts that cn+i = 1 — J2iLi c i an d J2iLi Ji = ~Jn+i, system (|7j) can be written as 
Vc' = — Aq(c') J'. The matrix Aq(c') defined in Section [2] is generally not symmetric and it 
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is not clear if it is positive definite. If Aj(c') is invertible (and we prove in Section [2] that 
this is the case), we can write system ([I])-© as 

(9) d t d - div(A (c')~ 1 Vc') = /(c) in Q, t > 0. 

Still, A (c') _1 may be not positive (semi-) definite. 

Our main idea to handle is to exploit its entropy structure. We associate to this 
system the entropy density 

N N 

(10) h(c') = 2Jcj (log ^ - 1) + cjv+i (log cjv+i - 1), ci,...,Cjv > 0, 2jcj < 1, 

i=l i=l 

where ctv +1 = 1 — J2iLi c « i s interpreted as a function of the other concentrations. Further- 
more, we define the entropy variables 

, . dh a 

(11) Wi = — = log-^, z = l,...,iV, 

and we denote by -ff(c') = V 2 /i(c') the Hessian of ft with respect to d. Then becomes 

(12) d t d - div(B(w)Vw) = /(c) in Q, t > 0, 

where u> = (wi, . . . ,wn) t and B{w) = A (c) _1 H (d)~ l is symmetric and positive definite 
(see Lemma ED- The advantage of the formulation in terms of the entropy variables is 
not only that the diffusion matrix B(w) is positive definite (which allows us to apply the 
Lax-Milgram lemma to a linearized version of (|T2|) ) but it yields also positive lower and 
upper bounds for the concentrations. Indeed, inverting (jTTI) . we find that 

(13) d = , i = l,...,N. 

v 1 1 + e Wl -\ he w " 

Therefore, if the functions W{ are bounded, the concentrations c« are positive and c * < 
1 which implies that c^+i = 1 — YliLi c« > 0. This observation is a key novelty of the 
paper. 

Formulation (jSJ) is needed to derive a priori estimates which are an important ingredient 
for the global existence proof. Differentiating the entropy H[c] = j n h(c)dx (now h(c) is 
interpreted as a function of all ci, . . . , cat +1 ) formally with respect to time, a computation 
(made rigorous in Lemma [H]) shows the entropy-dissipation inequality 

(14) ^ + K / V \V^\ 2 dx < 0, 

dt J n ^! 

where K > is a constant which depends only on (Dij). This estimate yields H 1 bounds 
for t/q. 

The existence proof is based on the construction of a problem which approximates f[T2"j) . 
We replace the time derivative by an implicit Euler discretization with time step r > and 
we add the fourth-order operator e(A 2 w + w), which guarantees the uniform coercivity of 
the elliptic system in V with respect to w. The existence of approximating weak solutions is 
shown by means of the Leray-Schauder fixed-point theorem. The discrete analogon of the 
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above entropy-dissipation estimate implies a priori bounds uniform in the approximation 
parameters r and e, which allows us to pass to the limit (r, e) — > 0. In particular, the 
entropy inequality provides global solutions. 

System ([I])-© with vanishing production rates, r = (r 1; . . . , r N+1 ) T = 0, admits the 
homogeneous steady state (c?, . . . ,(P N+1 ), where c° = meas(fi) -1 J n c®dx. We are able to 
prove that the solution, constructed in Theorem [TJ converges exponentially fast to this 
stationary state. For this, we introduce the relative entropy 

*+i r c . 
U*[c] = }] I Cilog^dx. 

Theorem 2 (Exponential decay). Let the assumptions of TheoremU\hold. We suppose that 
r = and minj =lj _jv+i ||c°||ii(n) > 0. Let (ci, . . . , c/v+i) be the weak solution constructed in 
Theorem[l\ and define c° = (dj, . . . , c° N+1 ) T . Then there exist constants C > 0, depending 
only on Q, and A > 0, depending only on Q and (D^), such that 

IIqO, t) - ^|Ui ( n) < Ce- x WH*[c% i = l,...,N + l. 

The proof of this result is based on the entropy-dissipation inequality (|14p and the 
logarithmic Sobolev inequality [16J, which links the entropy dissipation to the relative 
entropy, as well as on the Csiszar-Kullback inequality [26], which bounds the squared L 1 
norm in terms of the relative entropy. The difficulty of the proof is that the approximate 
solutions do not conserve the L 1 -norm because of the presence of the regularizing e-terms, 
and we need to derive appropriate bounds. 

The paper is organized as follows. In Section [21 we prove some properties of the diffusion 
matrices A{c) and A (c) and we show how system CO)-© of iV + 1 equations can be reduced 
to a system of N equations. Based on these properties, Theorems [1] and [2] are proved in 
Sections [3] and HJ respectively. For the convenience of the reader, the derivation of the 
Maxwell- Stefan relations ([2]) is sketched in Appendix |A] and some definitions and results 
from matrix theory needed in Section [2] are summarized in Appendix Q3J 

2. Properties of the diffusion matrices 

Let the Maxwell-Stefan diffusion matrix (D tj ) E ]R( Ar + 1 )x( Ar + 1 ) (jV > 2) be symmetric 
with Dij > for i ^ j and Da = for all % and set dij = l/D^ for % ^ j. Let c = (cj) G 
M. N+1 be a strictly positive vector satisfying YlfJi 1 c « = 1- We refer to Appendix [B] for the 
definitions and results from matrix analysis used in this section. According to (T2]) and (J7J), 
the matrix A — A(c) — (o^) G M( Af + 1 )x( Ar + 1 ) i s given by 

N+l 

aij = dijd for i, j = 1, . . . , iV + 1, j, an = - ^ dijCj for % = 1, . . . , N + 1. 

In p3l Section 7.7.1], the matrix with elements Qij / j Cj IS analyzed and it is shown that it 
is symmetric, positive semi-definite, irreducible, and a singular M-matrix as well as that a 
generalized inverse can be defined. Our approach is to apply the Perron- Frobenius theory 
to A, following [I]. 
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Lemma 3 (Properties of A). Let 5 = minj J=lv . .7V+1, i^j dij > and A = i^jdij- 
Then the spectrum cr(— A) of —A satisfies 

a(-A) c {0} U [5, A). 

The inclusion a (—A) C {0} U [5, 00) is shown in jH Section 5]. For the convenience of 
the reader and since some less known results from matrix analysis are needed, we present 
a full proof. 

Proof. The matrix A is quasi-positive and irreducible. Therefore, by Theorem [TBI of Perron- 
Frobenius (see Appendix [Bj , the spectral bound of A, s(A) = max{3ft(A) : A G cr(A)}, is a 
simple eigenvalue of A associated with a strictly positive eigenvector and s(A) > 9?(A) for 
all A e cr(A), A 7^ s(A). Here, 9R(z) denotes the real part of the complex number z. Thus, 

a {A) C {s{A)} U {z e C : < s(A)}. 

An elementary computation shows that c is a (strictly) positive eigenvector to the eigen- 
value A = of A. According to the Perron-Frobenius theory, only the eigenvector to s(A) 
is positive. This implies that s(A) = and 

a {A) C {0}U{z e C : M(z) < 0}. 

We can describe the spectrum of cr(A) in more detail. Let C 1 ^ 2 = diag( v /c7, . . . , ^/c^+i) 
be a diagonal matrix in $l( n+1 ) x ( n + 1 ) with inverse C _1//2 . Then we can introduce the 
symmetric matrix As = C~ x l 2 AC 1 ! 2 whose elements are given by 

a s f an if i = 1,...,N+ 1, 

\ d ijy /c~c~ = 1,...,N + 1, i^j. 

The matrix is real and symmetric since d^ = dji and therefore, it has only real eigen- 
values. Since A and As are similar, their spectra coincide: 

a(A s ) = a(A) C {0} U {z e R : z < 0} = (-00, 0]. 

Now, consider the matrix As (a) = As — a^fc ® y/c, where < a < 5 and sfc = 
(1/ci, • • • , y / c/v r +i) T . Then A s (a) is quasi-positive and irreducible (since a < 5 < cfy). 

Using Y^i=i °i = 1, a computation shows that —a is an eigenvalue of As (a) associated 
to the strictly positive eigenvector y'c. By Theorem [13] of Perron-Frobenius, the spectral 
bound of As(a) equals —a and 

a(A s (a)) c (-00, -a). 

Since As^a) and a^/c® y^c are symmetric, we can apply Theorem [T6l of Weyl: 

Xi(A s ) = X i {a^fc®^fc + A s (a)) < Xi{a^~c®^) + X N+1 (A s (a)), 

where % — 1, . . . , N + 1 and the eigenvalues Aj(-) are arranged in increasing order. Because 
of Ajv+i(As(cO) = —a and Aj(a;-y/c <g> -^/c) = for i = 1, . . . , N, Atv+i(«a/c <8> a/c) = a (see 
Proposition [14]), we find that Xi(As) < —a for i = 1, . . . , N and Xn+i(As) < 0. Thus, for 
all a < 5, 

a(A) = a(As) C {0} U ( 00, -a], 



A. JUNGEL AND I. V. STELZER 



implying that a{-A) C {0} U [5, oo). 

It remains to prove the upper bound of the spectrum. Denoting by 
norm, we find for the spectral radius of —A that 



\f the Frobenius 



-A)< 



A\ 



N+l 



'N+l 

\i,j=i 



1/2 




7V+1 



1/2 



d'ij 



A. 



<2 £ 

since < q < 1, finishing the proof. 



□ 



Lemma 4 (Properties of restrictions of A and As). Let A = A\i m (A) and As = As\i m (A s )- 
Then A and As are invertible on the images im(A) and im(v4,s), respectively, and 



(15) 



a{-A), a(-As) C [5, A), a((-A s )- L ) C (1/A, 1/5}. 



Proof. Direct inspection shows that ker{A) = span{c}, im(v4) = {1}^, where 1 
1) T G K.^ 4 " 1 , and ker(A s ) = spanji/c}. By the symmetry of As, it follows that 

(16) R N+1 = ker(A s ) ± © ker(A s ) = im(Aj) © ker(A s ) = im(A s ) © ker(A s ). 



Furthermore, using Theorem 
(17) 



since A = is a semisimple eigenvalue of A, 
R N+1 =im{A)®ker{A). 



We observe that both A and As are endomorphisms. Clearly, a (A) C cr{A) and a {As) C 
cr(As). We claim that is not contained in o~{A) or cr(v4s). Indeed, otherwise there 
exists x G im(yl) (or x G im{As)), x 7^ 0, such that Ax = (or Asx = 0). But this 
implies that x G ker(A) (or x G ker^s)) and because of ( fT7|) (or (TT6]) ) . it follows that 
x = 0, contradiction. Hence, A and are invertible on their respective domain, and ( TT51) 
follows. □ 

The above lemma shows that the flux-gradient relation ([7]) can be inverted since YlfJi 1 Ji 
= implies that each column of J is an element of {1}^ = im{A). In fact, we can write 
(13) as Vc = AJ and hence, J = A _1 Vc. Therefore, we can formulate (EL]) and (El) as 



d t c - div(-A _1 Vc 



r c 



in Q, t > 0. 



The next step is to reduce the Maxwell- Stefan system of iV + 1 components to a system 
of N components only. Still, we assume that Cj > for all i and X]^ 1 Q = 1. We define 
the matrices 



X 



l N+l 



/ \ 



W 



1 

Vo/ 



X 1 — Jjv+i + 





\1/ 



/A 
1 

Vo/ 



j(JV+l)x(AT+l) 
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where X -1 is the inverse of X and In+i is the unit matrix of ]R( Ar + 1 ) x ( Ar+1 ). A computation 
shows that 

where the (N x A r )-matrix Aq = (a° ) is defined by 



(19) a% 



n i — (dii - di,jv+i)ci iii^J, i,3 = !,■■■, N, 



Y,j=x,jfr(dij - di,N+i)cj + d itN+ i if i = j = 1, . . . , N, 



and the vector b = {pi) is given by bi = (i^jv+iQ) * — 1, . . . , N. In Lemma [5] below we show 
that A is invertible. Then, writing d = (ci, . . . , cat) t and J' = ( J 1; . . . , Jn) t , 

( V q c ^ = X - l Vc = (X- l AX)X- 1 J = (" A ° J ') • 

Thus, applying X~ l to ([IS]), since X _1 fltc = (<V, 0) T and X-V(c) = (r'(c),0) T , 

(20) <9 t c' - div(A l Vt/) = r' in ft, t > 0. 

We note that every solution c' to this problem defines a solution to (fT8j) and hence to 
©-(El) by multiplying (1201 (augmented via (d, 0) T ) by X and setting ctv+i = 1 — ^ i=1 Q. 

Lemma 5 (Properties of Aq). The matrix A e M ArxAr ; defined in (fl9l) . zs invertible with 
spectrum 

a(A )G[S,A). 

Furthermore, the elements of its inverse A^ 1 are uniformly bounded in Ci, . . . , £ [0, 1]. 

Proof. Since the blockwise upper triangular matrix — X -1 yLY is similar to —A, their spectra 
coincide and 

(21) a(A ) U {0} = <t{-X~ x AX) = a(-A) c {0} U [5, A). 

Observing that is a simple eigenvalue of —A, it follows that a(A ) C [5, A) and hence, 
Aq is invertible. 

It remains to show the uniform bound for the elements of Aq . By Cramer's rule, 

A^ 1 = adj(A )/ detA , where adj(A ) is the adjugate of Aq. The definition of Aq in ( TT9"j) 
implies that 

N 

\a%\ < ^ \ dik ~ d i, N+ i\ + \d i>N +i\ = Ki< K, i, j = 1, . . . , N, 

k=l,k^i 

where K = maXi = i j ... i jv^Q- Therefore, the elements of adj(v4o) are not larger than (N — 
l)!^ -1 . By (pTj) . the eigenvalues of A Q are bounded from below by 5. Consequently, 
since the determinant of a matrix equals the product of its eigenvalues, det(y4 ) > 8 . 
This shows that |a y | < (N - 1)\K N - 1 5~ N for all i, j. □ 
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Consider the Hessian V 2 h of the entropy density ffTUj) in the variables ci, . . . , c^. Then 

a) 



H = (hij) = V 2 h G R NxN is given by 



I 5iA 

"13 — 



hij = + — , i,j = l,...,N, 



where Sij denotes the Kronecker delta. The matrix H is symmetric and positive definite 
by Sylvester's criterion, since all principle minors det of H are positive: 

det H k = (ci • • • CkCN+iY 1 {^[1 Ci + ° N +^j > °> k = l,...,N. 

Lemma 6 (Properties of B). The matrix B = A^H' 1 is symmetric and positive definite. 
Furthermore, the elements of B are bounded uniformly mq,..., Cjy + i G [0, 1]. 

Proof. Using dij = dji and Yl!i=i °i = 1, a calculation shows that the elements of 
B- 1 = HA equal 



k,N+l . dik 

H c fc , 

Ci 



Pu = d hN+l (i- v Cfc Vl+J_)+ jr ( d -M± 

CiV+1 V fc=Wi / C7V+1 V *=Wi / *=Wij CiV+1 



where z, j = 1,...,N and i ^ j. Hence, 5 _1 is symmetric. We have proved above that 
H^ 1 is symmetric and positive definite. According to Theorem [151 the number of positive 
eigenvalues of A = H^B^ 1 equals that for B^ 1 . However, by ( |2T|) . A has only positive 
eigenvalues. Therefore, also B~ x has only positive eigenvalues. This shows that B~ x and 
consequently B are symmetric and positive definite. 

It remains to show the uniform boundedness of B. The inverse H~ l = (rjij) can be 
computed explicitly: 

_ f (1 - Ci)c,i Hi = j = l,...,N, 
lij \ -CiCj if i^j, i,j = 1,...,N. 

Denoting the elements of Aq 1 by ay, the elements fry of B equal 

N 

hi = oiu(l - Ci)ci - ^2 a ikCiC k , i = 1, • • • , N, 

k=l, kj^i 

bij = -auCiCj + aij(l - Cj)cj - ^ a ikCjC k , j, i, j = 1, . . . , N. 

k=l, ky^i,j 

By Lemma |5l the elements ay are uniformly bounded. Then, since q G [0, 1], the uniform 
bound for &y follows. □ 
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3. Proof of Theorem Q] 

The analysis in Section |2] shows that the Maxwell- Stefan system can be reduced to the 
problem 

(22) d t d - div(B(w)Vw) = r'(c) in Q, t > 0, 

(23) Vwi-u = ondQ, Wi(-, 0) = tw? in fi, i = l,...,N, 

where B = B{w) is symmetric and positive definite for ci, . . . , cjv+i > and q = Cj(u>) 
is given by ( fl"3l) . Furthermore, any solution c' = c'(w) to this problem defines formally a 
solution to the original problem ([I])-© by setting cn + i = 1 — J2iLi c i- We assume that 
there exists < rj < 1 such that c° > r\ for i — 1, . . . , AT and c^ +1 = 1 — $^ i=1 c° > 77. Then 
tu° = log(c°/c^ +1 ) satisfies u>° G L°°(fi), i = 1, . . . , N. 

Step 1: Existence of an approximate system. Let T > 0, m G N and set r = T/m, 
tk = rk for = 0, ... ,m. We prove the existence of weak solutions to the approximate 
system 

- / (c'(w k ) - c'(w k - 1 )) -vdx + [ Vv : B(w k )Vw k dx 
r Jn Jn 

(24) +e [ (Aw k ■ Av + w k ■ v)dx = [ r'(c(w k )) ■ vdx, veV N , 



where w k approximates w(-, tk) and e > 0. The notation ":" signifies summation over both 
matrix indices; in particular, 



/ 



N 

Vv : B(w k )Vw k dx = V / bijiw^Vvi ■ Vw k dx, 

A „-_1 J SI 



Q i,j=V 

and we recall that V = {u £ H 2 (Q) : V« • v = on (90}. The implicit Euler discretization 
of the time derivative makes the system elliptic which avoids problems related to the 
regularity in time. The additional £-term guarantees the coercivity of the elliptic system. 

Lemma 7. Let the assumptions of TheoremU\ hold and let w k ~ l G L°°(Q) N . Then there 
exists a weak solution w k G to ([21 



Proof. The idea of the proof is to apply the Leray-Schauder fixed-point theorem. Let 
w G L°°(n) N and a G [0, 1]. We wish to find w eV N such that 

(25) a(w, v) = F(v) for all v G V* 

where 

a(w,v) = / Vf : B(w) Vwdx + e / (Aw ■ Av + w ■ v)dx, 
Jn Jn 



a 



F(v) = J (c'(w) — c'(w x )) ■ vdx + a j r'(c(w)) ■ vdx. 
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Since B(w) is positive definite, by Lemma El the bilinear form a is coercive, 



a(w,w) = / 2_. bij(w)Vwi ■ Vwjdx + e I [\Aw\ 2 + \w\ 2 )dx > Ce\\w\\ 2 H 2, n -, N , 
■h'i j | Jn 

where C > is a constant. The inequality follows from elliptic regularity, using the 
assumption dQ G C 1 ' 1 . By Lemma [6] again, the elements of B(w) are bounded uniformly 
in c, and thus, a is continuous in V^xV^. Using < Ci(w), Cj(w fc_1 ) < 1 and the continuity 
of Ti, we can show that F is bounded in V N . Then the Lax-Milgram lemma provides the 
existence of a unique solution w G V N to ([25]). Since the space dimension is assumed to 
be at most three, the embedding H 2 (Q) >• L°°(f2) is continuous (and compact) such that 
w G L°°(Q) N . This shows that the fixed-point operator S : L°°(Q) N x [0,1] L°°(Q) N , 
S(w,a) = w, is well-defined. By construction, S(w, 0) = for all G L°°(fi) Ar . Standard 
arguments show that S is continuous and compact. It remains to prove a uniform bound 
for all fixed points of S{; a) in L°°(tt) N . 

Let w G L°°(p,) be such a fixed point. Then w solves (j25|) with u> replaced by w. 
Taking the test function v = w G V , it follows that 

- / (c(w) - c'iw^ 1 )) ■ wdx + / (Vw: B(w)Vw + e(\Aw\ 2 + \w\ 2 ))dx 
T Jn Jn 

(26) = a / r'(c(w)) • ifdx. 

in 

In order to estimate the first term on the left-hand side, we consider the entropy density 

h, defined in (TTOj) . Its Hessian is positive definite if ci, . . . , cn+i > and hence, h is convex, 

i. e. 

N N 

h{c) - h(c) < V/i(c) • (c - c) for all c, c G with < q, q, ^ c,-, ^ q < 1. 

3=1 3=1 

Using w = V/t(c'), we find that 

- / (c'(w) - c'iw"- 1 )) -wdx>- [ (h(c'(w)) - h{d{w k - x )))da:. 
T Jn T Jn 

By Lemma El B is positive definite: 

/ Vw : B(w)Vwdx > 0. 
in 

Finally, using the assumptions J2iLi r i( c ) = ~ r N+i(c) and J^fJi 1 r «( c ) l°g c i < 0> 

/ r'(c(w)) • wcfe = / n(c(w)) (log Cj(w) - logCAr+i(w)) dx 

Jn Jn \ i=1 J 

N+l 

rj(c(w)) logQ(w)(ix < 0. 



1=1 
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Therefore, ff26l) becomes 

cr / h(c (w))dx + er / (\Aw\ 2 + \w\ 2 )dx < o I /i(c(w; fc_1 ))rfx. 
Jn Jn Jn 

This yields an H 2 bound uniform in w and a (but depending on e and r). The embedding 
H 2 (Q) ^ L°°(Q) implies the desired uniform bound in L°°(Q), and the Leray-Schauder 
theorem gives a solution to ( |24|) . □ 

Note that we obtain the uniform bounds also under the weaker condition ([6]). In this 
case, fl26l) can be estimated as 



(27) a / /i(c'(w))rfx + er / (| Aw\ 2 + \w\ 2 )dx < a i h{d \w k ~ 1 ))dx + aC r rmeas{Q). 
Jn Jn Jn 

Step 2: Entropy dissipation. Since the diffusion matrix B(w k ) defines a self-adjoint 
endomorphism, the entropy-dissipation estimate 



Vw k : B(w k )Vw k dx > I X\Vw k \ 2 dx 



holds, where A is the smallest eigenvalue of B{w k ). Unfortunately, A depends on c(w k ) 
and we do not have a positive lower bound independent of w k . However, we are able to 
prove an entropy-dissipation inequality in the variables v/ 'ci(w k ) with a uniform positive 
lower bound. In the following, we employ the notation f(c k ) = (f(c k ), . . . , f(c% +1 )) T for 
arbitrary functions /. 

Lemma 8. Let w k G be a weak solution to fl24l). Then 



Vw k : B{w k )Vw k dx > i ( \vV^\ 2 dx, 

A Jn 

where c k = c{w k ) = (ci(w k ), . . . , CN + i(w k )) T is defined in ( TT3|) and C]y + i(w k ) = 1 — 

Proof. First, we claim that 

( Vw k : B(w k )Vw k dx = [ V\ogc k : (-A^V^dx, 
Jn Jn 

where A is defined in Lemma HI To prove this identity we set z 1 — (^1; • • • i %n) — 
B(w k )Ww k E JH Nxd and zjy + i = — J2i=i z i e ^- d - Then the definitions of w k and z^+i yield 

N 

Vw k : B(w k )Vw k = Vw k :z =Y^ (Vlogc, fc - Vlog4 +1 ) • Zi 

i=l 

N+l 

(28) =^Vlog Cl fc -^ = Vlogc fc :2, 

i=i 
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where z = (z',z N+1 ) T . Using Vw k = HVd(w k ) and B = A^H^ 1 , where H = H(c'(w k )) 
is the Hessian of h, it follows that z' = Aq Vcf(w k ) or, equivalently, Vc'{w k ) = A Q z'. A 
computation shows that for % — 1, . . . , N, 

N 

Vc k = (A z')i = ( d v - di,N+i)(ziC k - ZjC k ) + d itN+1 Zi = {-Az) t = (-Az) i} 

since each column of z is an element of im(A). Because of Az G im(v4), we have (—Az) jy+i = 
— J2i=i(~ Az)i = Vc k N+1 . We infer that Vc fc = — Az and consequently, z = (— A)~ l Vc k . 
Inserting this into (|28|) proves the claim. 

We recall from the proof of Lemma H] that the images of A = A\ im (A) and As = 
As\im(A s ) are given by im(v4) = {1} ± and im(A s ) = spanjv^}" 1 " = {C~ 1 ^ 2 x : x G 
im{A)}, where C+ 1 / 2 = diag((cf) ±1 / 2 , . . . , (4 +1 ) ±1 / 2 ) G R( N + 1 )*( N + 1 ) . Then the defini- 
tion -A = C^^-As)^ 1 / 2 implies that -A = C^i-Ag)^ 1 / 2 and hence, (-is)" 1 = 
C-^i-A)- 1 ^ 2 . We infer that 



Vlogc fc : (-I^Vc* = 4(VVc*) : C^-i^C 1 / 2 ^ 



c 



A- 



= iVVc 1 : (-i 5 ) _1 V\/? > — (Vv 7 ?! 2 . 

The inequality follows from Lemma H] since (—As) -1 is a self-adjoint endomorphism whose 
smallest eigenvalue is larger than 1/A. □ 

Step 3: A priori estimates. Next, we derive some estimates uniform in r, e, and r\ by 
means of the entropy-dissipation inequality. The following lemma is a consequence of (|26|) . 
the proof of Lemma [7J and Lemma [HJ 

Lemma 9 (Discrete entropy inequality) . Let w k G be a weak solution to (I24I) . Then 
for k > 1, 

H[c k } + ^- [ \V\^\ 2 dx + er [ (\Aw k \ 2 + \w k \ 2 )dx <U[c k -\ 
where c k = c(w k ) andH\c k \ = j Q h(c k )dx. Solving this estimate recursively, it follows that 

k k 

n[c k } + ^-J2 I \VVJ\ 2 dx + erJ2 [ {\^w J \ 2 +\w j \ 2 )dx<ri[c }. 



Let w k G be a weak solution to ( )24l) and set c k = c{w k ). We define the piecewise- 
constant-in-time functions w <yT \x,t) = w k (x) and c^ T '(x,t) = (c k ,...,c%) T (x) for x G f2, 
t G ((k — l)r, kr], k = 1, . . . , m, c^ T '(-, 0) = (c?, . . . , c%) T , and we introduce the discrete 
time derivative D T c^ = (c^ —a T c^ T ')/r with the shift operator (a T c^)(x, t) = c J \x 1 t—r) 
for x G SI, t G (r,T], (a T c (r) )(x, t) = c°(x) for i 6 (1, t 6 (0,r]. The functions (c (t) ,w (t) ) 
solve the following equation in the distributional sense: 

(29) D t c {t) - dw{A \c {T) )Vc {T) ) + e(AV T) + w {r) ) = r'(c (r) ), t > 0. 
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Lemma [9] implies the following a priori estimates. 

Lemma 10. There exists a constant C > independent of e, r, and rj such that 

(30) ||V / cW|| L 2(0,T;^l(n)) + \^\\w {T) \\ L 2^ T . H 2 m < C, 

(31) \\C {T) \\L'2(0,T;HHQ)) + || D T C {T) (I L2(0,T;V) < C. 

In the following, C > denotes a generic constant independent of e, r, and 77. 

Proof. Estimate ( 130]) is an immediate consequence of the entropy inequality of Lemma |9] 
and the L°°-bound for c^ T \ To prove (I3ip. we employ the Holder inequality: 

»l|2 _ A f T \\ A /jP)v7. fjr)\\ 2 



\^ C i \\l 2 (0,T;L 2 (Q)) - 4 / \\\/ C i ^\ C i \\L2(Q) dt 



T 

\\J n (r)\\ 

L°°(fi)ll v V ^ ll£ 2 (fi) L 



< J / ll\/cf T) llL^IIV\/ci T) || 2 r ,_tft 



< 4||cf ) || L oo (0iT;i oo (n)) ||Vycf ) ||^ 2(0)T . L2(Q)) < c, 

using ( )30l) and the fact that < < 1, z = 1, . . . , iV. Here and in the following, C > 
denotes a generic constant independent of s, r, and rj. By ( 1291) and L 2 (fi) ■=— >■ (if 1 (f2)) / , 



|^rC (T) || L2( o,T;V0 < C |K 1 Vc« || i2(0iT;L2 



(«)) 



+ £C||w (r) || i 2 (0iT . jf/2(f7)) + C 1 1 r ' (C (T) ) 1 1 ^2 ( ,T;i2 (fl)) 
< C||v4 1 ||z°°(0,T;L°°(n)) || Vc^ ||L2( ,T;La(n)) 

+ £:C||u; (r) || L 2 (0iT .^2 (o)) + C\\r'(c (T) )\\ L 2 {0iT . L 2 m . 

The proof of Lemma [5] shows that the elements of Aq 1 are bounded by a constant which 
depends only on iV and (-Df/). Since < < 1 and r' is continuous, (r'(c^)) is bounded 
in L 2 (0,T; L 2 (fi)). Therefore, in view of (JED and the bound on c (t) in L 2 (0,T; iT 1 ^)), 

II D T C^ || L2(0,T;V) < C7 

finishing the proof. □ 

Step 4- Limits e — > and r — > 0. We apply the compactness result of [12, Theorem 1] 
to the family (c^). Since the embedding H 1 ^) L P (Q) is compact for 1 < p < 6, (|3"T|) 
implies the existence of a subsequence, which is not relabeled, such that, as (e, r) — » 0, 



C W c' = (ci,...,cjv) strongly in L 2 (0,T;17(fi)), 1 < p < 6. 

As a consequence, q > 0, 5^=1 c « — 1' anc ^ c jv+i = 1 — Si=i c * — 0- Because of the 
uniform L°°-bounds for cf \ this convergence holds even in the space L q (VL x (0, T)) for all 
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1 < q < oo. Furthermore, by fl30]) - fl3Tj) . up to subsequences, 

Vc (t) Vc' weakly in L 2 (0, T; L 2 (fi)), 

D t c {t) -± d t c' weakly in L 2 (0, T; V) , 

e«; (T) ->■ strongly in L 2 (0, T; # 2 (ft)). 

Since the elements of Aq 1 are bounded and < c^' < 1, 

4)(c (t) ) -1 Ao(c')- 1 strongly in L 9 (0, T; L q (Q)) for all 1 < q < oo, 

r'(c (r) ) -» r'(c) strongly in L 2 (0, T; L 2 (fi)), 

setting c = (ci, . . . , ctv+i) t - The above convergence results are sufficient to pass to the 
limit (e, r) — > in the weak formulation of f[2Uj) . showing that c satisfies 

<V - div(A (c , )" 1 Vc') = r'(c) in L 2 (0, T; V). 

This proves the existence of a weak solution to (Q and ([3]) with initial data satisfying 
c® > i] > and 5^1=1 c° < 1 — 77. In view of the uniform bounds and the finiteness of the 
initial entropy, we can perform the limit 77 — > to obtain the existence result for general 
initial data with c° > and J2iLi c ? — 1- This proves Theorem [TJ 

4. Proof of Theorem [2] 

First, we prove that, if the production rates vanish, the L l norms of the semi-discrete 
molar concentrations are bounded. We assume that there exists < 77 < 1 such that c° > rj 
fori = l,...,JV + l. 

Lemma 11 (Bounded L 1 norms). Let r = 0. Then there exists a constant 70 > 0, 
only depending on c° , such that for all < 7 < min{l,7 } and sufficiently small e > 0, 
depending on 7, the semi- discrete concentrations c k = c{w k ), where w k G solves 
satisfy 

(32) (1 - 7 )l|c?||n(n) < ||c?|Ui(n) < (1 + 7)lk ? ll^(n), i = l,...,N, keN, 

N N 

(33) ||c^ +1 || L i(n) - 7^ ||c°|| L i(Q) < ||c^ +1 || L i(Q) < ||c^ +1 l| L i ( n) +7^ ||c°IUi(n)- 

i=l i=\ 

Furthermore, \\c k N+1 \\ L i {n) > ^\\c° N+1 \\ L i {n) > 0. 

Proof. We recall that r = T/m for T > and m G N. Let k G {1, . . . , m}. Using the test 
function v = e.{ in f )24|) . where is the i-th unit vector of ~R N , we find that 

/ c k dx = / c\~ x dx — st I w k dx, i — 1, . . . , N, 
Jn Jn Jn 

where we abbreviated c k = Ci{w k ). Solving these recursive equations, we obtain 
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Because of the £-terms, we do not have discrete mass conservation but we will derive 
uniform Z^-bounds. The entropy inequality in Lemma M shows that 

U[c k } +et [ ((Aw k ) 2 + {w k f)dx < U[c k } + er [ (\Aw k \ 2 + \w k \ 2 )dx 
Jn Jn 

(35) <H[c k -\ k>l. 

Solving this inequality recursively, we infer from H[c k ] > — meas(fi)(iV + 1) that 

k 

er Kll^(n) < ~ W] < W] + meas(fi)(iV + 1). 
i=i 

Consequently, using kr < T, 

1/2 



/ Wi\dx < erC V \\w{ < erCVk V 
j=i Jn i=i \i=i 

< C^eTk(K[c } + meas(ft)(iV + 1)) 



< Cy/eT(H[c } + meas(fi)(iV + 1)). 
Let 7 > and < e < 1 satisfy 

(36) 0< 7 <minM,7o= ( 2 V IK°|Un) I \\c\ 



N+l\\L' l (n) ( > 



(37) 0<VS< ^tk^Mm 

~ C^/T(n[c } + meas(n)(N + 1)) 



Then, in view of (J3 

(l-7)l|c?||i i i(n)<||cf|| i i( n ) = ||c?|| i i(n)-erX; / ^x<(l + 7 )||^ 

3=1 Jn 

These relations hold for all % = 1, . . . , N. For i — N + 1, we estimate (using (132 



on 



4 +1 rfx = jf M - g cA rfx > jf M - (1 + 7) g c° j 

„ v 1 

= / c^ +1 rfx- 7 V / c°cfe > -||c^ +1 || £ i (n) > 
Jn i=1 Jn A 



dx 



by definition of 70. A similar computation yields 



/ c k N+l dx< [ (l-(l- 7 )f;cnrfx=||c v+1 || Ll(n) +7f]||c°|| Ll (n). 
This proves the lemma. □ 
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For the proof of Theorem [21 we introduce the following notations: c k = (c k , . . . , c^ +1 ) T , 
where c k = meas(f2) _1 J n c k dx for i = 1, . . . , iV + 1, k > 0. Further- 



ed = c\ . . . ,c 



^k \T 



1, • • • , >-Ar + i; ? 



more, we set w k = (w±, . . . , w^) T , and u> fc = (w k , . . . , u^r) T , where w k = \og(c k / c k N+l ) for 
z = 1, . . . , N. We recall the definition of the relative entropy 



N+l 

H*[c k ) = £ / 
i=i Jn 

,k nT ,k : 



JV+1 „ fc 

c 

Ci log -^dx. 



Employing the test function u> — tu in ( liMl) . we obtain 

(c'(w fc ) - c'^^ 1 )) • - w k )dx + [ Vw k : B{w k )Vw k dx 

Jq 

+ e [ (\Aw k \ 2 + w k ■ (w k - w k ))dx = 0. 



We estimate the integrals term by term. 

Using the definition c k N+l = 1 — J2iLi c i-> a computation shows that 

(c'(w k ) - c^w*" 1 )) ■ w k = (c k - c*" 1 ) • logc fc . 

Therefore, we find that 

„k 

Jf„.M\ J/„..k-l\\ f„,.k -h i„ / ik 



(c'(w k ) - c'(w )) ■ (w k - w k )dx = (c k - c k ~ L ) ■ log —dx 
n Jn c 



[c k - c'- 1 ) • log ^dx + f{c k - c fc - x ) ■ log - k dx. 
n c° J n & 

The first integral on the right-hand side can be estimated by employing the convexity of 
h(c) as a function of Ci, . . . , cjv+i, which implies that 

h(c(w k )) - h(c(w k+1 )) < Vh(c(w k )) ■ (c(w k ) - cV^ 1 )) 

= log(c fc ) • {c k -c k - 1 ). 

Thus, because of J^i^ c \ = 1 an d ^ ne definition of H*[c fc ], 

/ (c k - c^ 1 ) • log ^dx > W[c k ] - H*[c k ~ 1 ]. 
Jn c 

For the second integral, we employ the bounds ( I32l) - (l33l) as well as 7 < 1 and e > 
sufficiently small, which yields 

1 1 A 7" 1 4 + l - 1 



< 4 < , i = l N, < < 

z n+i 



1 + 7 $ 1-7' ' 1 + 7 ci,, I-7' 
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where 7 = 7(l/c^ +1 — 1) > 0. Here, we have used again that J^^ 1 c° = 1. Then, with 
C\ = meas(f2), 

N „ 

7) da? 



f (c * - c*-i) . log £-dz > - / J> 4 fc log(l + 7 )dx- / 4 +1 log(l + 
c Jn i=1 Jn 

+ / J2 4' 1 log(l - 7)dx + / 4~\ log(l - 7)d 
Jn i=1 in 



>_ Cllog (i±ll(i±3). 
- ' s (1 _ 7)(1 __ } 

We have already proved that 

I Vw k : B{w k )Vw k dx >-r[ \VV^\ 2 dx. 
Jn A J q 

Applying Young's inequality to the e-term, it follows that 

H*[c k ]-n*[c k - 1 ] + ^ [ \vV^\ 2 dx<^ [ \w k \ 2 dx + C.log £ + 7 ffi + ?l 
^ Jn 2 Jn (1 -7JI 1 -7) 

The logarithmic Sobolev inequality [16J as well as the bounds ( 132|) -( 133|) show that 

N+ 1 r k N + l r 

W\c k \ = / c k log^dx+J2 / c k \og^dx 

i=l J V C i i=l J V °i 

N+l N 

<C(Q)J2 / \^Jc k \ 2 dx + Y, / cflog(l+7)<k+ / 4 +1 log(l + 7)dx 
i=l Jn i=1 Jn Jn 

<C(tt) [ |Vv / ?| 2 dx + C 1 log((l + 7 )(l + 7)), 
Jn 

from which we infer that 

(l + C 2 r)'H*[c k ]<'H*[c k ~ 1 ] + £ -^ [ \w k \ 2 dx + C^ 

* Jn 

where C 2 = 4/(C(0)A) and, for r < 1, 

c ^ = c ' log KiKi + c®a log((1 + 7)(1 + ^ 

We can estimate w by using the bounds for c k of Lemma [TTJ 

- jv - 

/ \w k \ 2 dx<J2 / (|logcf| + |log4 +1 |) 2 dx<C 3 , 
Jn i=1 Jn 

where C3 > depends on the L 1 -norm of c° and 7. Hence 

rT[c fe ] < (1 + Car)" 1 ^ [c*" 1 ] + yC 3 (l + ^r)- 1 + C 7 (l + C^r)- 1 . 
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Solving these recursive inequalities, we conclude that 

k k 

2 



n*[c k ] < (i + c 2 r)- k n*[c } + £ -^c 3 53(1 + cvr + c, 53(1 + c 2 r) 



3=1 3=1 
The sum contains the first terms of the geometric series: 

A- 



3 = 1 



yielding 



K*[c (T) (-,t)] < (i + c 2 T)-^n*[c°] + —1 + —L, t>o. 

2C 2 C 2 t 



Now, we choose sequences for e, r, and 7 such that 7 — >■ 0, C 1 /t — > 0, and (I37j) is satisfied 

(then also e — > 0). This is possible since C 7 — )■ as 7 — >■ 0. Then, because of c- r ^ — >■ q in 
L 2 (0, T; L 2 (n)) for z = 1, . . . , iV + 1, the limit (e, r, 7) -> leads to 

?T[c(-,t)] < e- C2 'H*[c°], t>0. 

Moreover, we can pass to the limit 77 — >• 0. Finally, since f n Ci(-,t)dx = f n c®dx for i = 
1, . . . ,N + 1 (see Lemma [TT1) . we can apply the Csiszar-Kullback inequality |26j to finish 
the proof. 

Appendix A. Derivation of the Maxwell- Stefan relations 

The Maxwell- Stefan relations (j5J) for an ideal gas mixture of A^+l components are derived 
by assuming that the thermodynamical driving force di of the i-th component balances the 
friction force We suppose constant temperature and pressure. Our derivation follows 
[1]. For details on the modeling, we refer to the monographs [151 127] - 

The driving force di is assumed to be proportional to the gradient of the chemical 
potential fa [271 Section 3.3]: 

(38) di =^Vfa, t = l,...,N + l. 

Here, R is the gas constant, T the (constant) temperature, and q = Pi/rrii is the molar 
concentration of the z-th species with the mass density pi and the molar mass m 8 of the 
z-th species. For more general expressions of di, we refer to [151 Chapter 7]. The chemical 
potential under isothermal, isobaric conditions is defined by fa = dG/dci, where G is the 
Gibbs free energy. In an ideal gas, we have Vfa = -RTVlogQ implying that di = Vq. 

The mutual friction force between the z-th and the j-th component is supposed to be 
proportional to the relative velocity and the amount of molar mass such that 

N+l , , 

(39) /* = - E CiCj[U ^~ Uj) , 

3=l>j¥>i %3 
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where D^ are the binary Maxwell- Stefan diffusivities [21, Formula (16)]. By the Onsager 
reciprocal relation, the diffusion matrix (D^) is symmetric [2U]. Then, using (|38p . ( IBTJj) . 
and the definition Jk = CkUk, the balance di = ft becomes 

N+1 r T -r T 

which equals (j2J). 

Another derivation of the Maxwell- Stefan equations (IT]) and (j2j) starts from the Boltz- 
mann transport equation for an isothermal ideal gas mixture [5J. The main assumptions 
are that a diffusion scaling is possible and that the scattering rates are independent of the 
microscopic velocities. Then, in the formal limit of vanishing mean-free paths, (JTJ) and (J2]) 
are derived. The diffusivities D^ are determined by 

_ T{mj + rrij) 

TTL j & 

where are the averaged scattering rates of the collision operator associated to the 
components i and j. Since Sij = Sji, the diffusion matrix (Dij) is symmetric. We also 
observe that (Dij) does not depend on the concentrations which is consistent with the 
assumption made in this paper. 



Appendix B. Matrix analysis 

We recall some results on the eigenvalues of special matrices such as symmetric, quasi- 
positive, or rank-one matrices. Although most of the results in this appendix are valid for 
matrices with complex elements, we consider the real case only and refer to the literature 
for the general situation [THJ [23], [25] . 

A vector x G R n (n G N) is called positive if all components are nonnegative and at 
least one component is positive. It is called strictly positive if all components are positive 
|25j . Let A = (ciij) G R nxn be a square matrix. The unit matrix in W axn is denoted by 
I n . Let cr(A) denote the spectrum of A. The spectral radius of A is defined by r(A) = 
max{|A| : A G a(A)}, and the spectral bound of A equals s(A) = max{9ft(A) : A G a(A)}. 
An eigenvalue of A is called semisimple if its algebraic and geometric multiplicities coincide 
and simple if its algebraic multiplicity (and hence also its geometric multiplicity) equals 
one. The following theorem is proved in [23l Theorem 3.4]. 

Theorem 12. Let A G W axn and let A G c(A) be a real eigenvalue. Then A is semisimple 
if and only if 

R n = im(A - XI n ) © ker(A - XI n ). 

The matrix A is called quasi-positive if A ^ and Ojj > for all i ^ j and irreducible 
if for any proper nonempty subset M C {1, . . . ,n} there exist i G M and j ^ M such 
that dji 7^ 0. If n = 1, A is called irreducible if A ^ 0. For quasi-positive and irreducible 
matrices, the following result holds [25j Theorem A. 45, Remark A. 46]. 
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Theorem 13 (Perron- Frobenius) . Let A be a quasi-positive and irreducible matrix. Then 
its spectral bound s(A) is a simple eigenvalue of A associated with a strictly positive eigen- 
vector and s(A) > 5f(A) for all A £ o~(A), A ^ s(A). All eigenvalues of A different from 
s(A) have no positive eigenvector. 

The spectrum of rank-one matrices can be determined explicitly [23| Section 3.8, Lemma 
2]. Notice that any rank-one matrix A £ ^ nxm can De written in the form A = x <E> y, 
where x eR n ,y e R m . 

Proposition 14 (Spectrum of rank-one matrix). Let x, y £ M. n . Then o~(x (g) y) — 
{0, . . . , 0, x ■ y}, i.e., x ■ y is a simple eigenvalue. 

We recall two results on eigenvalues of products and sums of symmetric matrices. 

Theorem 15. Let A £ W ixn be symmetric and positive definite and let B £ W axn be 
symmetric. Then the number of positive eigenvalues of AB equals that for B. 

For a proof, we refer to [231 Prop. 6.1]. 

Theorem 16 (Weyl). Let A, B £ R nxn be symmetric and let the eigenvalues Ai(A) of A 
and Xi(B) of B be arranged in increasing order. Then, for i — 1, . . . , n, 

\{A) + Xx(B) < \(A + B)< Xi(A) + X n (B). 
A proof is given in [18, Theorem 4.3.1]. 
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